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Abstract 

The production of charged particles has been measured in deep inelastic scattering (DIS) with the 
ZEUS detector. The evolution of the moments of the scaled momenta distributions in Q 2 and x has been 
investigated in the current fragmentation region of the Breit frame. The results in the current region are 
compared to e + e~ data and QCD analytical calculations. The results are consistent with the universality 
of single-particle spectra in DIS and e + e~ annihilation at high Q 2 . The results at low Q 2 disagree with 
analytical calculations based on the modified leading log approximation (MLLA) and local parton hadron 
duality (LPHD). 
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Introduction 



This paper reports the results of a study of the properties of the hadronic final state in positron-proton deep 
inelastic scattering (DIS). The event kinematics of DIS are determined by the negative square of the four- 
momentum transfer of the virtual exchanged boson, Q 2 = —q 2 , and the Bjorken scaling variable, x — Q 2 /2P-q, 
where P is the four-momentum of the proton. In the Quark Parton Model (QPM), the interacting quark from 
the proton carries four-momentum xP. The variable y, the fractional energy transfer to the proton in its rest 
frame, is related to x and Q 2 by y ~ Q 2 /xs, where s/s is the positron-proton centre of mass energy. 

A natural frame in which to study the dynamics of the hadronic final state in DIS is the Breit frame . In 
this frame the exchanged virtual boson is completely space-like and has a four-momentum q — (0,0, 0,— Q = 
—2xP Brelt ) = (E, p x , p y , p z ), where p Brelt is the momentum of the proton in the Breit frame. The particles 
produced in the interaction can be assigned to one of two regions: the current region if their z-momentum in 
the Breit frame is negative, and the target region if their z-momentum is positive. The advantage of this frame 
is that it gives a maximal separation of the incoming and outgoing partons in the QPM. In this model the 
maximum momentum a particle can have in the current region is Q/2. 

The current region in the Breit frame is analogous to a single hemisphere of e + e _ annihilation. In e + e~ — > qq 
annihilation the two quarks are produced with equal and opposite momenta, ±y/s/2. The fragmentation of 
these quarks can be compared with that of the quark struck from the proton which has outgoing momentum 
—Q/2 in the Breit frame. In the direction of this struck quark the scaled momentum spectra of the particles, 
expressed in terms of x p = 2p Brelt /Q, are expected to have a dependence on Q similar to that observed in 
e+e~ annihilation (|, ||[ || at energy *Js = Q. 

Within the modified leading log approximation (MLLA) there are predictions of how the higher order 
moments of the parton momentum spectra should evolve with energy scale J^, |(| . The parton level predictions 
depend on two free parameters, a running strong coupling, governed by a QCD scale A, and an energy cut-off, Qq, 
below which the parton evolution is truncated. The hypothesis of local parton hadron duality (LPHD) 0, which 
relates the observed hadron distributions to the calculated parton distributions via a constant of proportionality, 
is used in conjunction with the predictions of the MLLA allowing the calculation to be directly compared with 
data. 



Results 

The moments of the ln(l/x p ) distributions have been investigated up to the 4th order; the mean (I), width 
(w) , skewness (s) and kurtosis (k) were extracted from the distribution by fitting a distorted Gaussian of the 
following form: 

— - % ; oc exp (\k - \s8 - 1(2 + k)8 2 + ^sS 3 + ±-k8 4 ) 
a tot dln(l/x p ) *\8 2 4 V ' 6 24 j 

where 5 = (ln(l/x p ) - l)/w, over a range of 3 units (Q 2 < 160 GeV 2 ) or 4 (Q 2 > 160 GcV 2 ) units in m(l/x p ) 
around the mean. The equation was motivated by the expression used for the MLLA predictions of the spectra 
in ref. ||. 



Figure |l] shows the skewness of the ln(l/x p ) spectra as a function of ln(Q)- It is evident that the skewness 
decreases with increasing Q. Similar fits performed on e + e~ data shows a reasonable agreement with our results 
at high Q 2 , consistent with the universality of fragmentation for this distribution. The ARIADNE Monte Carlo 
model H gives a reasonable description of the data. The data are compared with the MLLA predictions of 
ref. [||, using a value of A = 175 MeV, for different values of Qo- The MLLA calculations predict a negative 
skewness which decreases towards zero with increasing Q in the case of the limiting spectra (Qq = A). This 
is contrary to the measurements. A reasonable description of the behaviour of the skewness with Q can be 
achieved for a truncated cascade (Qo > A), but a consistent description of the mean, width, skewness and 
kurtosis cannot be achieved || . A range of A values were investigated and none gave a good description of all 
the moments. We conclude that the MLLA predictions, assuming LPHD, do not describe the data. It should 
be noted though that a moments analysis has been performed [0 , taking into account the limitations of the 
massless assumptions of the MLLA predictions, where good agreement was found between the limiting case of 
the MLLA and e + e~ data over a large range of energy, 3.0 < y/s < 133.0 GeV. 
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Figure 1: Evolution of the skewness of the ln(l/cc p ) distribution with ln(Q). Data from e + e~ and ep are shown 
together with, ARIADNE Monte Carlo predictions and the MLLA predictions of Dokshitzer et al ^ predictions 
(the full line is Qo = A, the dashed Qo = 2 A, and the dotted Qq = 3A). The overlapping points are different 
x ranges in the same Q range. The inner error bars are the statistical error and the outer error bars are the 
systematic and statistical errors added in quadrature. 

Summary 

Charged particle distributions have been studied in the current region of the Breit frame in DIS. The moments 
of the ln(l/x p ) spectra in the current region at high Q 2 exhibit the same energy scale behaviour as that observed 
in e + e~ data. The moments cannot be described by the MLLA calculations together with LPHD. 
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